Transport in metallic multi-island Coulomb blockade systems: A systematic perturbative 

expansion in the junction transparency 
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We study electronic transport through metallic multi-island Coulomb-blockade systems. Based on a dia- 
grammatic real-time approach, we develop a computer algorithm that generates and calculates all transport 
contributions up to second order in the tunnel-coupling strengths for arbitrary multi-island systems. This com- 
prises sequential and cotunneling, as well as terms corresponding to a renormalization of charging energies and 
tunneling conductances. Multi-island cotunneling processes with energy transfer between different island are 
taken into account. To illustrate our approach we analyze the current through an island in Coulomb blockade, 
that is electrostatically coupled to a second island through which a large current is flowing. In this regime both 
cotunneling processes involving one island only as well as multi-island processes are important. The latter can 
be understood as photon-assisted sequential tunneling in the blockaded island, where the photons are provided 
by potential fluctuations due to sequential tunneling in the second island. We compare results of our approach 
to a P(E)-theory for photon-assisted tunneling in the weak coupling limit. 

PACS numbers: 73.23.Hk, 73.40.Gk, 72.70.+m, 85.35.Gv 
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I. INTRODUCTION 

Charging effects strongly influence the tunneling of elec- 
trons through small metallic, normal-state islands, leading 
to the well-known Coulomb blockade of transport at low 
temperatures^ 2 ^ Early theoretical and experimental studies 
explored Coulomb-blockade effects in a single gated island 
-a so-called single electron transistor (SET)-in great detail. 4 

Soon afterwards a number of theoretical 5 - 6,7,8 and exper- 
imental works also considered multi-island systems, where 
several islands are coupled capacitatively or by tunnel- 
ing junctions (see Fig. [TJ. One appro achSifliiiii^ is con- 
cerned with using one part of the multi-island system as 
measurement device for the residue. Serial arrays of 
islandsi2ii 4 -iiiiiJiiiilLi2i 2 i2i 2 i have been put forward for metro- 
logical purposes: for primary thermometry 16 and, oper- 
ated as electron pumps, as a standard for both current and 
capacitanceil 7 *!^!*^ For each cycle of periodically chang- 
ing the gate voltages one electron is transfered through the 
system. The main operation of these devices is described by 
first-order perturbation theory in the junction conductances, 
known as orthodox theory for single-electron tunneling. 22 
Contributions from second- or higher-order processes, so- 
called cotunnelingi^iSi are limiting the accuracy of the lat- 
ter devices.' 7 ^ In other experiments higher-order tunnel- 
ing contributions are of fundamental physical interest. For 
instance, logarithmical temperature dependence of the con- 
ductance associated with many-channel Kondo correlations, 
which govern the low-temperature transport properties of 
a single-electron transistori 2 ^ 7 ^ could be explained by a 
second-order perturbation expansion at the resonance peak. 29 
For large transparency of the tunnel junctions, first- and 
second-order perturbation theory is no longer applicable, and 
alternative approaches suitable to describe transport in the 
strong tunneling regime such as semiclassical approachesr^i 
real-time renormalization-group techniques, 8,32 path-integral 
Monte Carlo simulations, 33 or nonperturbative resummations 



within a real-time diagrammatic formulation 34 were put for- 
ward. 

In this paper we describe a real-time diagrammatic ap- 
proach to transport through multi-island systems that allows 
for a systematic perturbation expansion in the tunnel-coupling 
strengths of the tunnel contacts. In particular, we develop a 
computer algorithm that generates and calculates all possi- 
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FIG. 1 : An example of a multi-island geometry. The system consists 
of M = 5 islands with 1 = 4 island-island and J = 4 island-lead junc- 
tions. Any junction, e.g. i = 1, in between a pair (mj = l,m\ = 2) 
of islands is characterized by a transparency 0!q [see Eq. Jl H below 
for definition] and a capacitance C m= \ m '=2. Only some stray capac- 
itances are shown in the sketch. The state of the system is given by 
the excess charges on the islands: \%) = — 3, 1,2, —3,2) . 
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ble contributions to second-order transport for arbitrary multi- 
island systems. The theory is a generalization of a diagram- 
matic real-time approach that was invented for single-island 
devices^ In contrast to the latter, where a fully analytical 
treatment of first- plus second-order transport is feasible, 29 
the large number of second-order diagrams for multi-island 
systems motivates the development and use of a computer al- 
gorithm. The same idea has been used in advanced computer 
codes for SET networks based on orthodox theory. 35,36 

Second-order transport includes cotunneling processes in 
the Coulomb-blockade regime, where sequential tunneling is 
suppressed. In the standard description of these cotunneling 
processes* 2 ^ 4 , energy denominators appear that diverge when 
approaching the onset of sequential tunneling. These diver- 
gencies can be removed by replacing the energy denominator 
with some constant 25 or by partial resummation of higher- 
order contributions. 37 Besides cotunneling processes, there 
are other second-order contributions to transport that become 
relevant in the regime in which sequential tunneling is not 
suppressed. They account for the fact that quantum fluc- 
tuations due to tunneling give rise to a renormalization of 
both the charging energies and the tunnel coupling strengths. 
This results in transport contributions that have the functional 
form of sequential tunneling but with renormalized system 
parameters. An example is the tell-tale characteristics of a 
many-channel Kondo-effect at low temperature in a metallic 
single-electron transistor. 27 28 29 Similar renormalization ef- 
fects are also found in a diagrammatic real-time description of 
second-order transport through single-level quantum dots. 38 
The virtue of the real-time diagrammatic approach employed 
in this paper lies in the fact that the above-mentioned diver- 
gencies of energy denominators are automatically regularized, 
and that the renormalization effects are taken into account. It, 
therefore, allows for a complete and consistent evaluation of 
second-order transport. 

Second-order transport in multi-island systems is qualita- 
tively different from that in single-island devices. In multi- 
island systems, cotunneling processes in which two different 
islands change their charge occupation may occur. The energy 
conservation of the total process may be fulfilled by exciting 
one island on cost of the other one. This introduces a cou- 
pling between the two islands that may be of importance in 
devices in which one island is used as a measurement tool for 
the charge state of the other one, as we will discuss in more 
detail below. Our theory includes cotunneling involving either 
one or more than one island. 

The outline of the paper is as follows. First, we present 
the extension of the real-time diagrammatic theory 34 to multi- 
island systems. Then, in Sec.|lll| we present the algorithmic 
approach to generate and evaluate all second-order contribu- 
tions to transport. This includes a discussion of the applicabil- 
ity range of our perturbation expansion. Further details con- 
cerning calculation of diagrams and convergence properties of 
the perturbative expansion are included in AppendixeslAland 
IbI Afterwards, we illustrate our theory by applying it to one 
of the simplest multi-island setup, namely that of two single- 
electron transistors put in parallel, Sec. IIVI For this exam- 
ple, we discuss the physics of energy exchange between the 



two transistors, interpret it within a simplified picture using 
a P(E) -description of photon-assisted tunneling, and demon- 
strate the virtue and the limits of this picture by comparing it 
with our full second-order calculation. 



II. REAL-TIME DIAGRAMMATICS FOR METALLIC 
MULTI-ISLAND SYSTEMS 

A. Metallic multi-island geometries 

The system we consider in this paper consists of a num- 
ber M of small metallic islands as well as L leads, which are 
connected to some islands by J junctions (see Fig. Q. Ad- 
ditionally there are / junctions between the metallic islands. 
Accordingly with any junction i = 1, 2, . . . ,7 we will asso- 
ciate the pair of islands (;«,, mj) it is connecting, where the 
order is arbitrary but fixed, thus defining a 'direction' for the 
junction. Likewise for any of the junctions between leads and 
islands, 39 indexed by j = 1+ 1, I + 2, ...,/ + /, we define a 
pair ntj). Furthermore we have to take into account the 
capacitances C/ m between islands and leads, C mm i in between 
two islands, as well as capacitances C gm to additional external 
gates g. 

For typical samples of metallic islands the level spectrum 
is dense and the total charge on any island is large. Therefore 
not only the leads but also the islands can be considered as 
large equilibrium reservoirs, described by Fermi distribution 
functions. These equilibrium distribution functions are not in- 
fluenced by tunneling processes comprising only a very small 
fraction of the overall number of electrons. Electron-hole ex- 
citations left behind after tunneling are quickly equilibrated- 
on a time scale short compared to typical times between tun- 
neling events. 

The state of our system is then described by the excess 
charges % = {ni,---, «m} (total charge minus background 
charge) sitting on the M islands. 

The system can be modeled by the Hamiltonian: 

H = H l + Hm + V + H t =Ho+H t . (1) 

The noninteracting electrons in the leads and islands are 
described by 

L M 

H L = Y.H^v a J KV a l kv, Hm = Y* £ lv c l,Xv C mXv, 
1=1 kv m=\Xv 

(2) 

where wave vectors K and A numerate electron states within a 
given transverse channel v. (Note, that the transverse channel 
index v includes the spin of the electrons. For ease of notation 
we omit subindices jq v or A mv throughout.) 

Coulomb interaction of electrons is captured by V = 
V(x), the electrostatic energy of a given charge state \x) = 
\rii,ri2, ■ ■ ■ ,«m)- It depends in a complex way on gate and bias 
voltages and the resulting charges on the capacitances C; m , 
C mm r, and C gm , respectively. A straightforward scheme to cal- 
culate this dependence, governed by classical electrostatics, is 
given in Ref. |40j There a capacitance matrix is introduced 
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to first calculate voltages on the islands and subsequently the 
electrostatic energy. For an example of the gate and bias volt- 
age dependence of the electrostatic energy, see the discussion 
of a two-island setup in Sec. II V XI 

Finally charge transfer through the junc- 
tions is depicted by the tunneling Hamiltonian 



j kXv 



L L T !Vl ; ;ivS*'v e 

i XX'v 



I +C.C. 



(3) 



where and TV', , are tunneling matrix elements for junc- 
tions 7'and i, respectively, and exp(±/0 m// .) is a charge shift 

operator, acting on the charge state \%) described above^i 

Involving a phase operator </>„, as a canonical conjugate to 
the charge operator of the island m, i.e., [^ m ,N,n\ = i, these 
operators exp(±/0 m ) changes the excess particle number on 
the island m as n m — > n m ± 1 for each tunneling process ac- 
cordingly. 



B. Diagrammatic technique 

In this section we generalize the diagrammatic technique 
developed in Ref. [34] for a single SET to study multi-islands 
systems as described by Eqs. Q-O above. A short overview 
over the derivation is given, while the reader is referred to 
Refs. 29 34, and 42 for more details. In Sec. Ill B 51 we explic- 
itly discuss differences between single and multi-island case. 



1. The current 

One of the central objects of our theoretical description is 
the probability P x , to find the multi-island system in a certain 
state \%) = |«i,«2) • • • j«m) in charge space. Experimentally 
accessible quantities such as the average charge of a certain is- 
land follow directly from this probability. The other quantity, 
we will mainly be focused on, is the current flowing through 
junction j into reservoir given by the change in the number 
of particles, 



7,(0 = ej t (Nij(t)) 

= e L [TlM jKV c mjXv e- i ^J)(t)+ C . C .} 
kXv 



(4) 



Note that the current operator has a similar structure as the 
tunneling Hamiltonian Hj. 



2. Time evolution of operators 

The nonequilibrium time evolution of the charge degrees of 
freedom is described by the expectation value of the (diago- 



nal) density matrix 

Px(t) = ((\X)(x\)(t)) 



Tr (p rVV''^('')' \ X ) WT-e-W* 1 ^ 
Tr(po£H) s /*i [dt' 2 ... [dt' s 



t\ <tU<...<t' 



:T K [H T {t[) 1 H T (t l 2 ) I ...H T {t' s ) 1 \ X ){x\} . (5) 



Here we replace time-ordering T + and anti-time-ordering op- 
erators T by introducing integration / dt ' along the Keldysh 
contour with "times" t' running forward from fo to t and 
backward from t to fo (see Fig. [2j. The ordering of times 
t[ <t' 2 < . . . <t' s is with respect to this Keldysh contour with 
the Keldysh time-ordering operator T K arranging the opera- 
tors in the tunneling Hamiltonian in proper order. Hr(t)i de- 
notes the tunneling part of the Hamiltonian [Eq. Q] in inter- 
action representation with respect to Ho. 

Due to the separation of fermionic and charge degrees of 
freedom the Hamiltonian Hq, including interaction through 
the charging energy V{%), is bilinear only in electron oper- 
ators. Thus we can apply Wick's theorem and perform the 
trace over these degrees of freedom by contracting tunneling 
vertices in pairs. 

In fact, there are two independent contraction lines 
for each vertex, as each of the two electron oper- 
ators a]. KV /a ljKV and c f Xv /c m .x v constituting a tunnel- 
ing vertex is connected to another matching operator 
a ljKv/ a ]-Kv a - w ^ c m j kv/ c ] n xv' not necessar ily both of these be- 
ing part of the same tunneling vertex. However in the limit 
of large transverse channel number "simple loop" configu- 
rations dominate and we can represent contractions by one 
(directed 43 ) double tunneling line between pairs of vertices. 
Readers may refer to Figs. 2 and 10 of Ref.|34jfor an illustra- 
tion. 

A realization of a resulting diagram is shown in Fig. [5] 
where certain physical processes can be identified with parts 
of the diagram. The left-most part shows explicitly how the 
charge state is changed by a tunneling process across junction 
j. Processes with several tunneling lines at a time t correspond 
to higher-order processes such as cotunneling. 



3. Generalized master equation 

From these diagrams for the time evolution of the density 
matrix a formally exact master equation can be derived, 



dt 



Pz(*)= E / dt'[P x ,{t')I, x ,, x {t\t)-P x {t')Y, x , x ,{t',t)\, 



x'^x Jt ° 



(6) 

where the central object needed as input is the full quantum 
mechanical transition rate IL/ x (t',t) from a state %' at time t' 
to state x at time t. This rate is the sum over all irreducible 
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FIG. 2: Diagram describing the time evolution of the islands' 
density matrix with sequential tunneling, cotunneling (two tun- 
neling lines at a time), and higher order processes (from left to 
right). The diagram is reducible, i.e., it can be cut at a certain 
time without cutting a tunneling contraction line. Summing ir- 
reducible parts (shaded) yields the self-energy. In the left-most 
part a tunneling process through junction i in between a pair 
of islands (m,,mj) changes the state of the island on the up- 
per branch from \xo) = |«i,«2, •••,««!;, •••,««) to \xi) = 
\ni,ri2,..-,n mi - l,...,n m / + l,...,n M ) ■ 



( 1 ) Draw all topologically different diagrams with tunnel- 
ing lines and choose junction and direction for each line. As- 
sign energies V(x) to propagators, (0, to tunneling lines. 

(2) Each segment from t s < t < t s +\ gives a factor [AE S + 
if}] , where AE S is difference of left-going minus right-going 
energies (including the energies co r of tunneling lines). 

(3) Each tunneling line gives a rate function a- / '/' ± ((o r )- 
-as defined below-for the direction of the tunneling line 
across junction j/i going backward/forward with respect to 
the Keldysh contour. 

(4) There is a prefactor (—1) for each internal vertex on 
the backward propagator. 

(5) Integrate over the energies co r of tunneling lines. 

The golden-rule rate a ; (a !± correspondingly) stems from 
implicit integration over the energy of one of the double tun- 
neling lines 



(see Fig. diagrams and corresponds to the self-energy of a 
Dyson equation for the full propagator of the system. After 
Laplace transforming, the stationary distribution P x of charge 
states is found by 



= Yfx^x'^ with z x',x = [ \ dt> E z',z( f '>°)- 



(7) 



In a similar manner as for the density matrix, a diagram- 
matic representation for the expectation value of the current 
is found. The current operator becomes manifest in any such 
diagram as an additional external vertex [cf. Eqs. and 0] 
at the (right-most) time t . As above partial self-energies can 
be defined 



^C.0 = E{£g z (* , .0+^(* , ,0}- 



(8) 



a^(co) = j dE4f±(E+(0)fZ.{E) = ±<4 ± ^-ff_ ■ 

' (10) 

where f,, is the Fermi distribution function of lead / or is- 

J l/m 

land m and f~ = 1 — /+ . Hereby assuming constant tunneling 
matrix elements T ;v = as well as a constant temperature 
kgT = l/j8 across the sample and neglecting energy depen- 
dencies of the density of states N l / mv {E), each junction j is 
characterized by a single parameter, 



«0 = E |^ V |V' jV (0)A^ v (0) = 



1 



= Rk 

4n 2 e 2 Rj 4jc 2 Rj 



, (ID 



related to the tunneling resistance Rj of the junction. 



Here l/ , includes all those diagrams contributing to ^ x '. x 
with the rightmost vertex on the upper propagator and right- 
most contraction line describing tunneling out of lead lj 
through junction j as well as diagrams with right-most vertex 
on the lower propagator and the right-most process describing 
tunneling into lead lj ; correspondingly summing the remain- 
ing self-energy diagrams with right-most tunneling through 
junction j defines T, J , . 



The stationary current then follows as 



x,x 



x' L x'-x - ' 



(9) 



xx 



4. Diagrammatic rules 

Identification of terms of the sums in Eq. (|5jl with diagrams 
leads to a set of rules (for details see Ref. 34) for calculating 
the value of a certain realization of a diagram. Here we will 
give these rules in the form most convenient for calculating the 
Laplace transform of the self-energy, which allows calculation 
of stationary probabilities and current via Eqs. and (|9ji: 



5. Diagrams for single and multiple islands 

The formalism of real-time diagrammatics as presented in 
the previous section is applicable to the study of arbitrary 
multi-island system, as defined by the Hamiltonian Eqs. Q- 
(|3}. Hitherto derivation and application was mainly concerned 
with the case of a single island (see, e.g., Ref. l42l and refer- 
ences therein) and one study of two islands in series. 8 There- 
fore this section is devoted to a brief summation of differences 
between the single and multi-island cases. 

Virtually the entire formalism of real-time diagrammatics 
translates from the single to multi-island case by switching 
over from one variable of interest, the charge on the island n, 
to the M-tuple % = «i,«2, • • • ,«m of charges on all the M is- 
lands. This extension of state space has nontrivial effects only 
for the practical application of the method, but not for general 
scheme and derivation. Complications for the multi-island as 
compared to the single-island case arise due to the non-trivial 
electrostatic charging energies, numerous tunneling processes 
mediating between any two charge states, and most impor- 
tantly the exponentially increasing number of charge states as 
the number of islands grows. 
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III. ALGORITHMIC DIAGRAMMATICS 

As for a general multi-island system a closed analytic solu- 
tion seems unattainable, we set up and use in the following an 
automatized, computer-based numerical approach. 

A. Algorithmic scheme 

From our discussion above we learned about the essen- 
tial steps in our problem: Solving the electrostatics delivers 
the energy terms, needed for calculating rates of certain pro- 
cesses, which allow us to find stationary solutions of the rate 
equations. In this section we present the automatization of this 
process, where the (by far) most intricate part is the actual cal- 
culation of rates. This includes the automatic generation and 
calculation of real-time diagrams. 

The first step in our scheme is the restriction to a finite 
charge state space. This subspace will include all states within 
a few kgT around the classical ground state and addition- 
ally all states reachable by simple tunneling events from these 
classically occupied states. Note that this choice can be self- 
consistently checked from the resulting occupation probabili- 
ties. For all these states the electrostatic energies V{%), which 
also depend on gate and bias voltages, are calculated accord- 
ing to the method presented in Ref. 40. 



B. Transition rates 

We gain the rates entering Eqs. and up to second or- 
der in tunneling by generating and calculating all self-energies 
diagrams with one or two contraction lines. The scheme of 
generating all diagrams of a certain order is presented in the 
following for second-order diagrams, first order being trivial, 
whereas higher-order diagrams, although easily generated, are 
considerably harder to calculate and will not be considered be- 
low. 

Generating diagrams: Starting with a certain state %' of 
our chosen subspace of charge states for both, upper and lower 
branch of the Keldysh contour (see Ref.|2Jfor a discussion of 
the diagonality of E), the following steps were implemented. 

A first (i.e., left-most) vertex v\ is placed on either the up- 
per or lower branch, a junction ji/z'i as well as a direction 
for the tunneling event is chosen. This determines the charge 
state^i on both branches in between vertices vi and vi and 
consequently the energy terms determining the "free" propa- 
gation, governed by Hq and the tunneling line energy a>\ . For 
the second vertex V2 the same choices can be independently 
made, as it is not contracted with the first vertex. (Such a 
connection would result in a diagram, reducible to two first- 
order blocks, thus not being part of the irreducible second- 
order self-energy.) The third vertex, however, is connected to 
either one of the vertices V\ or Vi, inheriting whereby junction 
and direction of the tunneling event, whereas position on up- 
per or lower branch is free to choose. This is also the only 
freedom of choice left for the last vertex, which is combined 
with the remaining tunneling line. Just as at the first vertex, 



charge states have to be changed according to the tunneling 
processes at the vertices all along the propagator up to the 
right-most final state %, which is the same on upper and lower 
propagator. 

Estimating the number of diagrams allows some gauge of 
the complexity of the problem. For each charge state (for low 
temperatures three per island, resulting in a total of 3 M ) there 
are 2 4 x [2(J + I)] 2 x 2 second-order diagrams. Here the first 
factor stands for the upper/lower branch for each vertex, the 
second results from vertices Vi and \>2 for choosing junction 
and direction of tunneling, while the last factor stems from 
combining vertices into pairs. 

For each of these diagrams an analytical integral expression 
is immediately given by the diagrammatic rules, discussed 
above. Evaluation of these expressions to a numerical value is 
discussed in AppendixlAl 



C. Stationary state solutions 

The value of any first- or second-order diagram is then 
added to the appropriate matrix element of the self-energy ma- 
trix EyYy > where indices %' ,% are initial and final charge state 
of the diagram. 

All diagrams required for calculating currents are already 
created within this scheme by identifying the right-most tun- 
neling vertex with an external current operator. Accordingly 
adding (with the proper sign) the values of diagrams, where 

the last vertex involves junction j, yields £yfy • 

x! x cons ti mte me fr^t terms of an expansion of the self- 
energy in powers of Oq, 



E , - V E l 

k=\ 



x'x 



(12) 



(k) 

Expanding the probabilities P x — YX=Q^x correspondingly, 
we find from the stationary rate equation in first and second 
order, 



x' x' x' 



which finally gives solutions for P£/ and P^,, For first- and 
second-order current 45 , we likewise find 



4" = -"L^H'S 1 

x,x' 



(14) 



,(2) 

/} = -ie 



x,x 



Thus we derived a scheme to calculate experimental accessi- 
ble quantities for an arbitrary multi-island geometry, namely 
the average current in any lead. 
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D. Limits of the approach 

There are both practical and fundamental limits of the 
approach presented above. As the number of charge 
states, and consequently the number of diagrams to calcu- 
late grows rapidly with the number of islands, a straight- 
forward simulation of some interesting existing experimen- 
tal applications 17,46 using long arrays of islands (7-100 is- 
lands) is practically infeasible. Similarly, including higher 
than second-order contributions would involve more complex 
integral expressions that complicates the numerics consider- 
ably. 

The fundamental limit deals with the question of conver- 
gence of the perturbative expansion. As we derive in Ap- 
pendix IB1 the inequality max{£#r, |A — fij } > |cr(A)| has to 
be fulfilled for tunneling across any junction j. Here |A — fij\ 
is the distance to resonance for a tunneling process across 
junction j with chemical potential difference flj and d(A) 
is the self-energy, characterizing the spectral density of the 
level into which tunneling occurs at energy A. This incor- 
porates renormalization of level position as well as a finite 
lifetime width, as reflected in real and imaginary part of the 
self-energy. We elaborate on this statement in Appendix 151 
where we analyze the example of a single SET by comparison 
to a nonperturbative approach. 



IV. AN EXAMPLE: TWO PARALLEL 
SINGLE-ELECTRON TRANSISTORS 



To give an example we will now apply our method to calcu- 
late transport through two single-electron transistors (SETs), 
where the two islands are capacitively coupled to each other, 
see Fig. |3] This system was studied experimentally, e.g., in 
Refs. Il0| and |47| The capacitive interaction implies that a 
change of the charge state of one SET will change the effec- 
tive gate voltage of the other, which will lead among other 
things to a broadening of the charge state transitions as de- 
scribed in Ref. [nj This is an effect visible already in first 
order of the tunneling conductances. The second-order con- 
tributions include cotunneling in each SET separately, but also 
second-order processes in which both SETs are involved. For 
the latter, the energy gained in a tunneling process in one SET 
can be used to excite the other SET to a charge state not ac- 
cessible for first-order transitions. This results effectively in 
an energy exchange between the two islands, a feature that 
is qualitatively new as compared to the single-island case or 
first-order transport in multi-island systems. We concentrate 
our analysis on a regime, where it is this energy exchange, that 
enables transport in an otherwise blockaded region. One may 
try to simulate this energy transfer in a simplified picture of 
photon-assisted first-order tunneling in the second SET cou- 
pled to a energy-providing photon bath that represents elec- 
trostatic fluctuations by tunneling processes in the first SET. 
This simplified approach and the effect of the energy-transfer 
processes on the transport characteristics will be described in 
more detail below. The advantage of this approach is its fea- 
sibility and compact analytical results. The disadvantage is 
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FIG. 3: Two capacitively coupled SETs. The upper one is the noise 
generator, and the lower the quantum noise detector. 



that it relies on several approximations, and its applicability 
range is unclear. This question can be answered, though, by 
our full-fledged second-order transport calculation. By com- 
paring the results of both approaches we are able to define the 
range of parameters beyond which the simplified picture loses 
its reliability. 



A. Setup 

The system consists of two metallic islands (called noise 
generator and detector as explained below), each connected 
to a separate source and drain lead, see Fig. [5] Each junc- 
tion is generally described by a capacitance and a tunnel re- 
sistance, as discussed in Sec.|ll] For simplicity we here as- 
sume that for each island the connection to source and drain 
is symmetric, and also that the source-drain voltage is applied 
symmetrically. Note, that we change our notation compared to 
the general setup in Sec. Ill Al in order to make the association 
with the noise generating and detecting part more transpar- 
ent. We denote the generator/detector junction capacitances 
Cq/Cd, the resistances Rq/Ro, and the applied source-drain 
voltages V G A /V^ d . Furthermore the working point of the gen- 
erator/detector SET is controlled by gate voltages (Vq/Vq) ap- 
plied across the gate capacitances (Cg/Cf,). The two islands 
are coupled through the coupling capacitance C c . 

Following Refs.|4(Jand[lO we find the electrostatic energy 
of our system. This is most straightforwardly described intro- 
ducing the sum capacitances of the two islands 



-G/D — 2C G/D 

the rescaled charging energies 



C 8 

^G/D 



■Cc, 



* 2{Cl-Cl/Cl) 
the interaction energy 



£int — 



e 2 C c 



2{<% -<%/(%)' 



r L r z - r2 ' 



and also the dimensionless gate charges n^ D = C G ^ D V^ D /e. 
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Then the electrostatic part of the Hamiltonian reads 



V{n Gl n D ) 



E G (n G -n G x f+E D (n D -n D x f 

+£int("G-"*)("D-«?)- (16) 

B. The noise detector SET 



To illustrate the effect of photon-assisted tunneling clearly 
we consider our setup biased so that sequential tunneling is 
blocked in one of the SETs (the noise detector). This is done 
by gating n x < 1/2 such that the Coulomb energy difference 
A = Ed (I — 2n® ) between the ground state no = and first 
excited state «d = 1 is larger than temperature and applied bias 
A > {eV D d /2,k B T}. Neglecting the coupling to the other SET, 
the direct current through the detector is given by cotunneling 



eV D d 



1 



h 24 7T 2 \R D 



(eV«Y 



{2nk B T) 



A 2 



(17) 



where we disregarded cotunneling through the state nr> = — 1. 
We will find that the current induced by photon-assisted tun- 
neling can be made orders of magnitude larger than I® ot . 



C. The noise generator SET 

The other SET (the noise generator) is biased such that a 
substantial current is produced by sequential tunneling, but 
for simplicity we keep the bias low enough so that only two 
charge states are involved. To be specific n x 1/2, and 
eV G d /2 < 2Eq. The current through the noise generator SET 
is 



eV^R 



K 



seq 



4h R G 



„-eVg/2k B T 



-1 



(18) 



Each electron tunneling in the forward direction will gain the 
energy eV G d /2 from the applied bias. Generally the electron 
will dissipate its excess energy on the island or in the reser- 
voirs, but in our setup there is also the possibility to give the 
energy to a tunnel event in the detector SET. 



D. The effect of interaction 

The low-frequency effect of the interaction between the 
SETs can be described as an effective gate charge determined 
by the charge state of the other SET 

n x,e«( n D) = n° - (n D - n D x ) , 

and 



as discussed in, e.g., Ref . fiol and fl2l We want to minimize 
this effect, in order to clearly show photon-assisted tunnel- 
ing. That means choosing the detector bias low enough so 
that sequential tunneling is exponentially suppressed for both 
effective gate charges «f eff (0) and «f eff (l) ■ 

We will now estimate analytically the detector current due 
to photon-assisted sequential tunneling, driven by photons 
emitted from the generator with an energy higher than the 
Coulomb gap of the detector. We can then compare these ana- 
lytical estimates with the numerical results from our algorith- 
mic diagrammatics. We use P(E) -theory, 48 considering the 
noise generator SET as the environment of the noise detector 
SET. 

The function 



P(e) 



dtexp (j(t)- 



e 



(19) 



expresses the probability of exchanging the energy £ with a 
certain environment in a single tunnel event. Here the function 



7(f) = ([0(f)-0(O)] 0(0)} 



tr 



1 



CO" 



(20) 



is given by the correlator of the phase fluctuations on the de- 
tector island induced by the electron number fluctuations on 
the generator island. To simplify the analysis we approximate 
the fluctuations to be Gaussian. Then their asymmetric noise 
spectral density is given by 



(CO) 



dte- ia "(8N G (t)8N G (0)). 



(21) 



We want to calculate the photon-assisted rate for an electron 
to tunnel onto the detector island 



01 



IRk 
hR D 
\_Rk 
hR D 



= -— / de 



dede'f(e) [l -f(e' - A±)] P(e - e' 
eP(-e-A±) 



I _ e -£/k B T 



(22) 



where A± = A ± eV D d /2, and A± is the effective energy gap for 
an electron tunneling onto the island from the right/left lead. 
Thus we need P(e) for e < — A± indicating the probability to 
absorb an energy larger than A± from the environment (see 
Fig. 01. In this regime, noting that A± ^S> Ej nt and using the 
short time expansion e J ^ = 1 +J(t) we can approximate 49 



P(e) 



^int Sn( £ ) 

h e 2 • 



e <C -Ei, 



(23) 



In the relevant frequency regime and at zero temperature the 
generator SET noise spectrum is^S^LS 



J N 



(*) = * 



r(e) 



7i4r(0) 2 



(24) 



where 



Q 

n °eff (»g) = - 7^ ("G - n° ) 

L G 



m = -— 

2n R G 



eV, 



sd 



£ < 



(25) 
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cited state is 



FIG. 4: Sketch of functions relevant to calculation of the photon- 
assisted tunneling rate [Eq. <22H . Without energy exchange with 
the environment [P{e) is 8 peaked, dotted line] the photon-assisted 
rate is given by the rate function a~ (broken line) at the effective 
energy gap A±. Coupling to the noise generator SET changes P(s) 
(solid line) and the rate is given by the integration in Eq. 1261 . The 
inset shows the asymmetric charge noise of an SET as calculated in 
Ref.lH. 



Furthermore S$(e) = for £ < — eV G d /2, which is related to 
the fact that one cannot extract more energy from a single tun- 
nel event than what is given by the bias voltage. 

We thus find the photon-assisted rates for electrons to 
tunnel onto the detector island across the right/left junction 
(cf. Fig.l), 



1 01 



1 Rk 
2%h R D 



eVg/2-A ± 



de 



er(-e-A d 



(£ + A±) 2 4r(0) 2 



:+A±)2- 
(26) 

Here the lower limit of the integral comes from putting tem- 
perature to zero in Eq. \22l and the upper limit from the fact 
that P(e) = for £ < -eV G d /2. Making the generator bias 
substantially larger than the detector gap (eV G d ^> 2A±), we 
find the approximate photon-assisted rates 



where 



,p±_ eVg RkRkEI 
01 4%K 2 hR D R G A| 



/(2T(0)/A ± ), 



(27) 



f(x) = ^ Mi - 2arctan - +xln (1 +xr) - 2x 



(28) 



We note that these rates are proportional to the bias voltage 
of the noise generator. If the detector gap is not too small we 
have r(0) <C A± which implies i«l and f(x) approaches 
the limit /(0) = 1 . For simplicity we assume this limit from 
here on. The relaxation rates of the detector are 53 



D± _ R K A ± 
10 ~Rd h 



(29) 



and accordingly the probability to find the detector in its ex- 



P\ } = 



i 01 



Rk 



ev G c int 



TD+ 
1 10 

R K eV. 



T D- 
1 10 

G c int 



R c 487z; 2 (A + +A_) 



1 

Ai 



1 

A 2 " 



R G 487T 2 A 3 



(30) 



and the photon-assisted tunneling (PAT) detector current (as- 
suming eVp -C A±) 



R K R K eVgeVg% 



■DrD-\ 



J A 1 : 



1 01 



■2 
'int 



R d Rg 327T 2 /iA 3 



(31) 



Comparing this with the standard cotunneling current in the 
detector, given in Eqn. d!7i . 



I D 

1 cot 



3 R D eVg ( Ei 



4R G A 



e ysd 



(32) 



we find that the photon-assisted detector current can be made 
substantially larger than the usual cotunneling current by us- 
ing high resistance tunnel junctions in the detector and a weak 
detector bias. 



E. Results of diagrammatic technique 

With the diagrammatic technique developed above, we 
have a method at our disposal for analyzing the complete two- 
island system on an equal footing. We do not perform a sep- 
aration into detector and noisy environment as with the P{E)- 
theory. In particular, we can study the mutual influence of the 
two SETs also for the strong coupling case. 

In the following, we will nonetheless focus on the weak 
coupling case; here we can clearly show the effects of PAT 
and sensibly compare to the results of P(E)-theoiy [Eq. OIH . 
Low-frequency effects, as described by the effective gate 
charges, and back-action are all intrinsically included in the 
calculation, but parameters are chosen to minimize these. 

Figure [5] (solid lines) shows current in the detector SET 
upon changing the bias on the generator SET for four differ- 
ent sets of sample parameters and gating. The detector is gated 
close to degeneracy point, i.e., «^ eff («d = 0) = 0.5, as the de- 
tector is set preferentially unoccupied in the Coulomb block- 
ade region. The detector is to remain in this blockade regime, 
independent of the generator's charge state. To achieve con- 
siderable, experimentally detectable, current through the de- 
tector, we must not be "too deep" in Coulomb blockade. 
Hence, we move from a region of well-defined blockade in 
curve A of Fig.|5]to a "mixed" regime in curve D (£;„( ps A±) 
by changing gating (and coupling). The exact parameters are 
given in the captions of Fig. [5] 

For small driving of the generator, we find the standard co- 
tunneling current in the detector SET. The cotunneling cur- 
rent shown in Fig. [5] (dotted horizontal lines for curves A-C) 
is calculated for a single SET with effective gating s {hq), 
weighted with the probabilities for the two generator states 
J'(ii g = 0)w1/2wP(iig = 1). 
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FIG. 5: Results of the diagrammatic technique for the genera- 
tor/detector setup (solid lines) compared to standard cotunneling 
(dotted lines) and photon-assisted current with slope as predicted 
from P(E) theory (dashed lines, x-intercepts fitted). All curves at 
T = 25 mK, eVjf = 0.5 k B K, and n° ef[ {n D = 0) = 0.5 with a J Q = 0.01 
for each junction; Eq = 10£#K, Ed = 10fc#K for A-C and 5fc#K for 
D, E- m = 0.5k B K for A-C and 1 k B K for D. as follows: A — 0.3, 
B — 0.4, C — 0.45, D — 0.4. 



For strong driving the condition for PAT above is ful- 
filled (i.e., eV G d > 2A±) and we find the detector current in- 
creasing proportional to the driving voltage of the generator 
[cf. Eq. J27H . For the curves A and B we find the slope as 
predicted by P(£')-theory [Eq. OH . dashed lines in A-C]. 
For curves C and D we do not meet all assumptions made in 
the derivation of the PAT rates. Consequently PAT and low- 
frequency interaction effects intermingle in a complex man- 
ner and the result deviates from Eq. 13 H . For D the cou- 
pling is chosen so strong, that back-action effects are also rel- 
evant, i.e., a meaningful definition of effective energy gaps 
and, therefore, comparison to PAT and standard cotunneling 
results is not possible anymore. 

As the photon-assisted rate was calculated assuming strong 
generator driving only, the crossover between standard cotun- 
neling and PAT is not properly described. The results of our 
full theory (solid lines in Fig.[5J show, however, the expected 
behavior-an onset of PAT (i.e., the x-intercept of the dashed 
lines in Fig. Incurves A and B), where eV^ ~ A± and a cor- 
responding crossover, when the PAT rate reaches the standard 
cotunneling rate. 

In conclusion our theory for the generator/detector setup 
shows standard cotunneling for low generator bias and 
photon-assisted tunneling for high generator bias, where the 
crossover depends on the effective energy gap of the detector. 
For weak generator-detector coupling our results agree per- 
fectly with a P(E) calculation of photon-assisted sequential 
tunneling, treating the generator as an energy providing envi- 
ronment of the detector. While the relative effect is strong- 
PAT current is of the same order as the standard cotunneling 
current-the overall currents are quite small. Stronger coupling 
yields higher currents, but now current is not given by photon- 
assisted sequential tunneling alone and our simplified P(E) 
calculation is insufficient. 

In the example of a detector-generator setup above we 
choose parameters to a regime where we can isolate one par- 



ticular new type of cotunneling process (involving energy ex- 
change between islands) as dominant constituent of current. 
Then we could compare to a P(E) calculation of photon- 
assisted sequential tunneling, simple enough to achieve an- 
alytical and physical transparent results. To reach such simple 
analytical results we made use of a number of assumptions 
and approximations: Within the P{E) calculation we did not 
consider standard cotunneling through the detector island (nor 
through the generator island). In calculating the noise spec- 
trum of the generator determining P(E) we neglect its cou- 
pling to the detector, which is permissible for appropriately 
chosen parameters. The Gaussian approximation was used 
connecting the autocorrelation of phase fluctuations on the 
detector island to the asymmetric charge noise of the gener- 
ator, while higher moments were neglected. The influence of 
higher moments (and consequently their detection) has been 
studied in a number of other systems recently. 54 ' 55 ' 5 ^2i^iSTo 
achieve analytically tractable results we exercised a short time 
expansion, which is justified in the weak coupling regime. Fi- 
nally, the zero temperature and large bias limit was assumed. 

The latter two approximations are in no way crucial and 
for convenience and simplicity of results only. Reaching 
such transparent results for comparison to the perturbative ap- 
proach is the main purpose of the P(E) analysis we performed 
here. Relinquishing this notion of achieving simple analytical 
results, a P(E) analysis may also be further refined to loosen 
some of the restrictions imposed above: Consider cotunnel- 
ing through the detector island, which depends on the state of 
the generator island (by the effective gating), which is occu- 
pied with either zero or one electron with some probability. 
This results in a shift of the Coulomb gap of the detector is- 
land dependent on the generator state. A P(£)-type theory 
of cotunneling in the detector can be set up to incorporate 
that effect of the environment. This will then just correspond 
to the "standard cotunneling" results in Fig. [5] which takes 
this very effect into account. In principle a P(E) analysis can 
also be employed without making use of the short time expan- 
sion to extend towards strong coupling. However, such an ap- 
proach is limited by the increasing importance of backaction 
effects in the strong coupling regime. Indeed, by construc- 
tion backaction effects are not included in a P(£')-theory. Its 
basic principle is to describe the environment by a single func- 
tion [namely P(E)], which in consequence cannot depend on 
the system's state dynamics. Backaction effects are, however, 
fully accounted for within the perturbative approach, where 
there is no separation into system and environment parts as in 
P(E). 

It should be noted that we implemented a perturbative ap- 
proach up to second order in the coupling. In the scenario 
above, the main contribution to current stems from sequen- 
tial tunneling in the detector (proportional to the detector 
coupling) assisted by fluctuations due to sequential tunneling 
through the generator (proportional to the generator coupling). 
By this way of counting powers of the coupling, it is obvious, 
that a second-order analysis cannot capture the influence on 
detector current of higher-order quantum fluctuations in the 
generator. 

We have demonstrated that cotunneling contributions in- 
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volving tunneling in different islands, where energy is ex- 
changed between electrostatically coupled SETs, can signif- 
icantly contribute to transport. A normal state SET is there- 
fore inherently sensitive to finite frequency noise of charge, 
to which it couples electrostatically. A number of other sys- 
tems have recently been suggested and used as on-chip mea- 
surement devices for finite frequency current and/or voltage 
noise: a single 60 or double quantum dot, 49 a superconducting 
qubit, 61 or a tunnel junction, both in the normal and super- 
conducting s tateM£2£££2i2££2& A normal state SET is in- 
herently sensitive to integrated noise of a broad energy win- 
dow, due to the peculiarities of the sequential tunneling rates 
[see Eqs. i22\ and J26H . In contrast, other systems mentioned 
above are particularly designed to make use of sharp reso- 
nancelike features in the detector to have high sensitivity for 
noise detection of a certain frequency. 



(a) / dujiai(±uji) J rfw 2 a2(±^2) 



W2, 



x 1 1 1 

(b) J <ko\aii±LOi) J e?cj2a 2 (±w 2 ) 
-„ -„ 

,'' ^2 ,'' X 



1 



FIG. 6: Two topologically distinct types of diagrams are found. This 
yields four different integral terms, to be evaluated. 



V. SUMMARY AND CONCLUSIONS 



In this paper we extended the real-time diagrammatics of 
Ref.Qto model transport through metallic multi-island sys- 
tems. We discussed our approach to automatically generate 
and compute diagrams in executing a systematic perturbation 
expansion up to second order in the tunneling conductance. 
This corresponds to sequential and cotunneling terms. Con- 
vergence properties of the perturbative expansion were ana- 
lyzed. 

In a setup of two coupled SETs, we demonstrated the im- 
portance of cotunneling involving both islands, where energy 
can be exchanged in between the two tunneling processes. 
This is linked to the notion of photon-assisted tunneling. We 
performed a P(E) analysis, treating one SET as an environ- 
ment for the other one and found excellent agreement between 
both approaches in the relevant weak coupling limit. In this 
regime one SET works as a detector of the finite frequency 
charge noise of the other SET. 



1. Reducing denominators 



The general structure of integral expressions for a diagram 
follows easily from the diagrammatic rules discussed above. 
A self-energy diagram of order R features R integrals over 
contraction line energies co r with corresponding rate functions 
a j± (co r ) but also (2R - 1) denominator factors [AE + /J]] -1 
for each segment in between two vertices. The workhorse for 
evaluating these expressions is Cauchy's formula, 



1 lex 

— — =P--inS{x), 
x+ir\ x 



(Al) 



which obviously can only be used when integrating over the 
variable x. Consequently we have to reduce the number of de- 
nominator terms from (2R — 1 ) to R. Note here, that Cauchy's 
formula can still be applied to higher-order poles, 
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APPENDIX A: CALCULATING DIAGRAMS 

In this appendix we explain how integral expressions, cor- 
responding to the diagrams created, are efficiently calculated. 
First of all, the number of expressions we have to evaluate can 
be reduced by identifying several distinct types, only differ- 
ing in some parameters. For these small number of terms nu- 
merical and analytical methods of automated calculation are 
discussed. 



While this reduction is possible for arbitrary order, we will 
illustrate it for second-order diagrams only. As illustrated in 
Fig- EI there are two topologically distinct diagrams, depend- 
ing on whether the first vertex is connected to third or fourth 
vertex. Note, that the position of vertices on the upper or lower 
branch is irrelevant for the general structure of the analytical 
expressions, which are indicated in Fig. [6] Elementary al- 
gebraic manipulation of the denominator factors-merely the 
treatment of infinitesimal rj requires some care-yields four 
expressions with the number of denominator factors reduced 
as required. 
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From the expression in Fig.|6ja) we get 
1 1 1 



©i — Ai + if] a>i + © 2 - A 2 + irj a>i - A 3 + irj 

' A A d \ 1 

Ai = A 3 



= < 
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Correspondingly, Fig.|6jb) yields 
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1 
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1 \ 1 



©1 + ©2 - A 2 + if] ) A[ + A 3 - A 2 
1 1 



©1— Ai + z'Tj © 2 -A3 + /?7 
for the second topologically distinct type. 

2. Mirror rule 



Further simplifications arise from an inherent symmetry 
property of the diagrammatic rules. From the construction 
principles for diagrams it is easily seen that, starting from any 
diagram, we will obtain another possible diagram by the fol- 
lowing operation: 

Reflect the diagram along a horizontal line, whereby ex- 
changing the forward and backward branches of the contour, 
and change the direction of all tunneling lines. The charge 
state for any part of the contour then remains unchanged. 
Both diagrams contribute to the same matrix entry of the self- 
energy, as only the diagonal part of £ has to be considered. 
Application of the diagrammatic rules shows that merely the 
energies in the denominator terms change sign, 



1 



1 



AE + iri 



-AE + ir\ 



1 



AE + iri 



(A3) 



resulting in cancellation of the real parts of the two mirrored 
diagrams. Consequently we only calculate the imaginary part 
of each diagram. This allows to utilize the 5-function part 
of Cauchy's formula for evaluating one of the two integra- 
tions. Our problem reduces then to the calculation of a single 
(principal value) integral, where the resulting integrand can be 



written as a product of one or two rate functions-there might 
be a derivative acting on one of them-multiplied with the prin- 
cipal value term P~. 



3. Evaluating the integrals 

Above [see Eq. ( I10H we calculated analytical expressions 
for the rate functions for the case of normal reservoirs with a 
constant density of states (DOS) and equal temperatures for 
both sides of the tunneling barrier. It is interesting, though, 
that important general features of the rate functions will per- 
sist, if these conditions are relaxed, e.g., for nonequal temper- 
atures, relevant for considering thermopower or self-heating 
effects, or for tunneling rates of quasiparticles in supercon- 
ducting devices, where the gapped BCS density of states en- 
ters. Common to all these cases, however, is the existence of a 
certain onset energy, where the asymptotic behavior is linear 
on one side, and vanishing on the other side of this thresh- 
old; the asymptotic convergence is hereby governed by some 
Boltzmann factor. 

Evidently, a high-energy cutoff is needed for convergence 
of integrals, where the two rate functions grow in the same 
direction. This cutoff is provided for in a natural way by the 
bandwidth of the reservoirs, as a more careful analysis of the 
rate function for a DOS limited to some finite bandwidth re- 
veals. The cutoff will at most enter logarithmically into any 
final result-indeed, for the single SET, it is known to drop out 
completely from the measurable quantities current and aver- 
age charge. Hence, there is no need for a microscopically 
detailed derivation of a specific cutoff function and we can 
choose for convenience of calculation a Breit-Wigner factor 
centered around the onset. 

An analytical solution of the integrals is then possible for 
the rate functions as given in Eq. dim . Contour integration 
leads to sums over Matsubara frequencies, which in turn re- 
sults in analytical expressions involving Digamma functions, 
where the onset energies enter as parameters (cf. Ref. l29l) . 

Since further algebraic manipulation of these complex 
terms is not feasible for multi-island geometries, it can be 
more convenient to adopt a numerical approach, already for 
calculating the integrals-the additional benefit being that this 
approach is also capable of effortlessly dealing with the alter- 
native rate functions discussed above. Numerical evaluation 
is helped by the general features shared by all these rate func- 
tions. They allow for precise numerical integration for a small 
region (of a few kgT) around the onset, while the asymptotic s 
can be trivially integrated analytically. 



APPENDIX B: CONVERGENCE OF A PERTURBATIVE 
EXPANSION 

In this appendix we address the question of the applicability 
range of our perturbation expansion. To motivate a criterion 
for general multi-island systems let us first consider the case 
of a single SET. In this case, we have a nonperturbative result 



12 



for the current within the so-called resonant tunneling approx- 
imation (RTA)2i at hand, 
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where a r = a+ + and the propagator 11(a)) = 
[ft) — A — (7(a))] -1 dressed with first-order transitions as ex- 
emplified in the self-energy 
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where the cutoff energy Eq is of the order of the charging 
energy. We see that the dominant contribution to the current 
comes from energies around the charging energy A. Higher- 
order quantum fluctuations embodied by <j(co) both shift (real 
part) and broaden (imaginary part) the region of contributing 
energies. 

For definiteness we consider now the scenario fi^ — jJ.R = 
eV > 0, where we are interested in the onset of sequential tun- 
neling around iil < A, while a finite bias eV 3> ksT is applied 
asymmetrically across the SET. The main contribution to the 
current in this regime is then given by electrons slowly tun- 
neling onto the island from the left lead (the rate for this bot- 
tleneck process being a^(co), while they quickly tunnel off to 
the right lead [o^(ft)) ^ (ft))]. Consequently the current 
[Eq. dBl H simplifies to 
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where the imaginary part of the self-energy, the lifetime width 
of the resonant level, is dominated by the (decay) rate for tun- 
neling out to the right lead. 

To understand the effect of a systematic perturbation ex- 
pansion of the current operator, as performed in the main part 
of the paper, we expand now this expression for the current 
[Eq. JB4H in orders of the tunneling conductance Oo. Expan- 
sion of the denominator yields 
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where the denominator terms are regularized by infinitesi- 
mal imaginary parts. This should be construed in terms of 
the Cauchy identity for generalized functions and derivatives 
thereof [Eqs. (|ATJ and iA"2tl. 

The term of this expansion of order a ( " contains the de- 
nominator (a) — A + /rj)~(" +1 ', which results in a n th deriva- 
tive of a 5-function. Within the integral expression for the 
current [Eq. JB4H we will use partial integration, so that the 
derivatives act on (co)o(co)" . The expansion in terms of 



|o(A)| 



)( fails 




A CO ^ 



FIG. 7: Sketch of a^(co) (broken line) and ImYl(co) (solid line) 
in the integrand of Eq. <B4L Perturbative expansion works, if the 
rate a£ is approximately linear in the region given by the peak in 
Imll(ffl), corresponding to |<r(A)| < ksT on resonance (a) or to 
|A — ft£, \ > |<t(A)[ in (c) and (d). The expansion fails (b) for ImII(o) 
incorporating the kink in a£. 



Oq therefore correspond to a Taylor expansion around the bare 
resonance A. In Fig.0we visualize the integration of Eq. JB4> 
for different parameter regimes. We integrate the rate func- 
tion a + (a>), which has a kink at CO = fit with a characteristic 
width of kgT, multiplied with Imn(ffl), which is peaked close 
to A, shifted by | Re a ( CO ) \ and has a width given by | Im a ( CO ) \ . 

As just shown, the perturbation expansion consists of ex- 
panding the peak of Imn(o)) in terms of derivatives of 5- 
functions, consequently capturing properies of 05 + (ft)) (and 
its derivatives) at A only. 

As is intuitively clear from Figs.0a)-0d), this will work, 
if the peak Imn(o)) is sharp, as compared to the structure in 
a^(co) in the relevant regime. This means that either |cr(A)| 
has to be smaller than ksT for considering the system on res- 
onance A = fit [Fig. 0a)] or that the expansion takes place 
so far away from resonance with the left lead, that the peak 
Imn(ffl) is not incorporating the kinked region of 0C + (ft)); 
i.e., |A-ju L | > |ct(A)| [see Fig. 0c) and0d)]. Then the 
rate function is expanded in the asymptotic region, where it is 
exponentially suppressed [Fig. 0c)] or rising approximately 
linearly 0d)]. Consequently higher-order derivatives from 
higher-order terms of the expansion in ao do not contribute 
and second-order perturbation theory suffices. 

We can formalize these intuitive arguments by considering 
the Taylor expansion of the bottleneck rate, resulting from 
the derivatives of the 5-function in the perturbative expansion 
[Eq. ffia i. 



CIqCO 



exp(j3o)) - 1 



CO 

2" 



jSco 2 j8 : 



CO"* 



6-2! 30-4! 



Note the increasing powers of the inverse temperature j3 , en- 
dangering convergence for low temperatures. The largest term 
now to appear in the n th term of the ao expansion is a term 
with magnitude 



d" m - l a+(co)a(co) 



n-l 



n-\ 
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where all derivatives have acted on 0^(0). 63 Correspond- 
ingly we get a series in the parameter /3 |(j(A)|, which has 
to be small for convergence, confirming our argument above. 

In summary, we found a criterion for the applicability range 
of second-order perturbation theory. The ideas developed by 
considering a single SET can be easily generalized to an ar- 
bitrary multi-island geometry. For each island, all tunneling 
rates should be approximately linear in the region of the con- 
tributing energies. This leads to the condition 

max.{k B T,\A-fij\} > |ct(A)|, (B6) 

where d(A) is the self energy characterizing shift and width of 
the resonance due to quantum fluctuations. For kgT > d(A) 
the tunneling rate is approximately linear on the relevant scale 



of integration, while for |A — /Xy > (7 (A) the kink in the rate 
function is outside the integration region, and consequently 
the rate is either linear or exponentially small. 

An example of the failure of perturbative expansion is an 
SET at finite bias and low temperature. Likewise in a setup of 
two islands in series 21 the resonance is shifted by a self-energy 
scaling as Re a °= (XqEc, not vanishing for aligned levels of 
the two islands. This can result in a failure of perturbative 
expansion at low temperature s£ 

No problems are encountered for the scenario laid out in 
Sec. II VI Here the generator SET is completely dominated by 
sequential tunneling, while the detector, biased on the order of 
temperature, correspondingly yields a small |cx(A)| sa a^ksT, 
fulfilling Eq. dBrfl . 
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